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1 Introduction 

In many practical problems, it is important to construct the set of all possible 
initial states of a dynamical system for which the trajectories do not violate 
given phase restrictions. For example, one of the basic problems in acceleration 
technology concerns maximal capture of particles. It is equivalent to maximiza- 
tion of the domain of stability [Q. On the other hand, the practical stability of 
motion generalizes the well-known Chetayev's definition of {A, A, to, Tj-stability 
Q . Analysis of stability on a finite time interval is realized using Lyapunov func- 
tions method. In works have been proved the sufficient conditions 
for practical stability and, for example, in |l|] , || , JtJ , Q have been obtained the 
necessary conditions in particular cases. The similar results hold for parametric 
and discrete systems [Q , [|| • 

The purpose of this paper is twofold. First, we shall prove the necessary 
and sufficient conditions for practical stability of nonlinear dynamical system at 
general phase restrictions. It will be shown, that in such a case the Lyapunov 
function is nondifferentiable. But if the set of initial data is starry compact, 
then it is possible building Lyapunov function which belongs to differentiable 
functions class. Our second goal is to estimate the optimal sets of initial condi- 
tions in structural forms for linear system and concrete phase restrictions. This 
estimations are important in theoretical sense and enable numerical methods to 
be devised [Ml. 



2 Auxiliary propositions 

Let p be a metric in the space R n . We shall denote by dK the set K frontier 
and by intK the interior of K, K C BP. The following proposition is true. 
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Lemma 2.1 If K C R n is a compact then there exists a function a : R n —* R , 
a G C (R n ), such that: a(x) = 1 at x G dK ; a (x) < 1 at x G intK ; a (x) > 1 
at x & R n \intK. 

Proof. As it is known, the function ip (x) = min p(x,y) is continuous. 

Moreover ip (x) > at x G i? n and ?/> (x) = if and only if x belongs to dK. Con- 
sider the function a\ (x) = 1 — ip (x) , x G if. This function is continuous on K, 
ai (x) < 1 if x G mf-K" and ai (x) = 1 if x G On the other hand the function 
a 2 (x) = l+ip (x) belongs to C (R n \intK), a 2 (x) > 1, x G R n \K and a 2 (#) = 1 

if x G Then the function a (x) = I a } J if > x satis- 



a 2 (x) = 1 + V (x) , a; G i?™\A: 
ties the lemma's conditions. 

Denote ||x|| the Euclidean norm of vector x G R n , S r (a) = {x G R n : 
||x - a|| = r}, A r (a) = {x G A" : ||x - a|| < r}, L x — {u ^ R n : u — Ax, A > 
0}, x G i?™\{0}. 

Definition 2.1 TTie set A C _R" is said to be starry if it is closed, G intE and 
for arbitrary x G R n \ {0} either L x <Z E or there exists a unique point z G L x 
so that z G BE. 



b G A then it is possible to build the curve z (A) = I , Q .' \ ^ ?' ' , and 



The starry sets are singly connected. In fact, if A is a starry set, a G A and 

Aa, A G [0, 1] , 
-Ab,A G [-1,0] 
z (A) g A, A g [-1, 1], z (1) = a,z (-1) = 6. 

If A is a starry compact then: a) for any I G Si (0) there exists a unique 
a > such that a/ G 9A and a = max k ; b) Minkowski's function fc (A, x) = 

max fc is continuous and positive, x G iZ n \ {0}. Really, the set E is 



bounded. Therefore the proposition a) follows from definition BTl Let us prove 
b). Consider x„, x G i? n \{0), lim x v = x, L = irnr, Z = Trfrr- There exists 

p^oo ll^pll 1^1 

a unique point z p G 9 E so that -p 2 ^ = l p . Inasmuch as d E is compact, then 
we can choose a subsequence {z m } C lim z m = u and u G <9A. Hence 

rn — >oo 

lim lp = I = jrV. Taking into consideration that k (E, x m ) = ||z m ||, k (E, x) = 

p—>oo II'" II 

||u|| and ||z TO || — ► ||tt||, m — > oo we have k(E,x m ) — > k{E,x), m — > oo. This 
proves b). 



3 Necessary and sufficient condition for prac- 
tical stability 

Let us consider a system of differential equations 



dx 

Tt =/(*,*) (1) 
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where a; is a n-dimensional vector, the vector function / (x, t) satisfies the con- 
ditions of the existence and uniqueness theorem, / (0, t) = 0, t € [to, T]. 

We denote x(t) — x(t,Xo,to) the trajectory of the system (|l|) at Cauchy 
condition x (to) = xq and also assume that <3?t C R n are compacts, G $i, 
te[t ,T], G C$ to , 0eG o . 

Definition 3.1 The unperturbed solution x (t) = of the system (j?p is said to 
be {Go,&t,to,T}-stable if x (t, xq, to) G $ t , t G [to,T] as soon as xo G Go- 



Theorem 3.1 For the trivial solution of the system fijy to be {Go, $t, to, T}- 
stable it is necessary and sufficient that there exists a continuous nonincreasing 
on the system (ji^ solutions Lyapunov function V [x, t) such that 

{xeR n :V(x,t)<l}C$ t ,te[t ,T}, (2) 

G C {x G R n : V(x,t ) < 1} . (3) 

Proof. Necessity. Let us G* be the optimal by inclusion set for {Go, $t, to, T}- 
stability of the unperturbed solution x (t) = 0. It means that zero solution of the 
system ([[]) is {G», $ t , t , T}-stable and if the solution x (t) = is {D , t , T}- 
stable then Dq C G*. Further, the set G* is compact. Following by lemma 
[2.l| one can build a function a G C (R n ) such that a(:c) = 1 at x G 9G», 
a(x) < 1 at i G intG* and a (x) > 1 at x G R n \intG*. Let us consider 
the function V (x,t) = a (tp (t, t , x)) where tp(t,to,x) — Xq on the solution 
x (t, Xq, to). Since V (x (t, xq, to) ,t) — a (xo) the function V (x, t) is nonincreas- 
ing on the system (|l|) solutions. We shall prove (|2|) by contradiction. Let there 
are x G R n and t G [to,T] for which the condition (|^) is true, but x ^ $ t . 
Then V (x,t) = (a (tp (t, to, x))) = a (x ) < 1 and x G G». Therefore 
x = x(t,xo,to) G <£>t. We obtain a contradiction. The validity of (H) follows 
from relations G C G* = {a; e R n : a (x) < 1} = Lr e R n : V (x, t ) < 1}. 

Sufficiency. Suppose that conditions (||) and (||) are satisfied but there are 
r G [to,T] and xo G $t such that the trajectory x(t) = x(t,xo,to) leaves 
the set $ T . From (^) follows inequality V (x (r) , r) > 1 which contradicts the 
condition (|^) . Our assumption is false and the assertion of the theorem is true. 

Corollary 3.1 If the trivial solution of the system is {Go,§t,to,T}-stable 
then the function 

{1 + min p (ip (t, t ,x) ,y) ,x G {x G R n : tp (t, t , x) G i?"\G*} , 
yedG, 
1- min p(tp(t,t ,x) ,y) ,x G {x G R n : tp (t,t ,x) G G*} 



satisfies all the theorem 3.1 conditions. Here p is a metric equivalent to Eu- 



clidean one. If Gq is compact then it is possible to take the Lyapunov function 
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in the form 



{1 + min p (cp (t, t , x) , y) , x E {x E R n : ip (t, t ,x) E R n \G } , 
vedG 
1 - min p (t, t , x) , y) , x E {x E R n : ip (t, t , x) E G } ■ 
y£dG 

4 Starry sets of initial conditions: criterion for 
practical stability of dynamical system 

We assume the set Go C R n is a starry compact. Let us introduce a new variable 

J fe(Go.xo) ' x ° + °> and denote yo = { k(Go,xo)> Xo ? °' if XQ u o then 
[ 0,a;o = [ 0,a;o=0. 

= f(k(G ,y )y,t)k- 1 (G ,y ). 

T3 ... / . n f f (k(G ,y )y,t)k- 1 (G ,yo) ,y ^ 0, . 
Putting g(y,t,y ) = 4 n n we obtain 

^ U, yo — U 

the system of differential equations 

^ = 5(y,*,M)).tG[to,r] (4) 

with initial condition y (to) = yo € K\ (0). 

Thus, the problem of {Go, $t, to, Testability of the system (|l]) zero solution 
can be formulated as the problem of {K\ (0) , to, Testability of the system 
(g) trivial solution, where * t = {(y,y ) ■ yk(G Q ,yo) & $t}, Vo G i^i (0). From 



theorem 3.1 follows the next 



Proposition 4.1 Tor t/ie trivial solution of the system fy) to be {Ki (0) , ^t,to, T}- 
stable it is necessary and sufficient that there exists a Lyapunov function V (y, t, yo) 
continuous and increasing on the system (Qj solutions and 

{(y, y ) : V {y, t, y )<l}Q^t,tE [t , T] , 

Ki{0) C {y :V(yo,t ,yo) < 1} ■ 

If p is Euclidean metric then the Lyapunov function can be taken in the form 

V (v> t, yo) = min s p ( z , <p (*: *o, y, yo)) • 

zGSi(O) 

Here cp (t, to, y, yo) is the system (Q) first integral in Cauchy's form. 
Let us consider the function 

fx Xo 
W(x,t,xo) = V[ — -,t, 



k (Go, a;o) ' k (Go,x ) 
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mm p[ z, <p t,t , — — r, — — 

zeSi(o) V V k{G 0l x ) k(Go,x ) 



■2-0 



V Mo, 



fc(G ,a;o) /c(G ,xo) 



Theorem 4.1 The trivial solution of the system is {Go,$? t ,to,T}-stable if 
and only if there exists a continuous Lyapunov function W (x,t, xq) such that 

{(x, x ) ■ W (x, t, x ) < 1} C {(x, x ):x£<Pt,x £Go},t£ [to, T] , 
G <={x :W (x , t 0: x ) < 1} . 



5 Practical stability estimations for linear sys- 
tem 

We assume that the set of initial conditions is a ball, that is Go = K c (zq), and 
phase restrictions are the sets 

T t = {xeR n : \lj(t) x\ < l,s= 1,2,... ,N},t £ [t ,T]. 

Here l s (t) is a continuous n-dimensional vector function, s = 1,2, . . . , N , 
parameter c is positive, zq £ R n . 

Following by corollary 3.1 Lyapunov function can be chosen in the form 
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- „ mm „ \\v{t,to,x) - y\\ ,x e {x : ip{t,t ,x) G R n \G } 
y-\\y-zo\\=c 

1- min \\(p(t,t ,x) - y\\ ,x e {x : (p(t,t ,x) e G } 



V(x,t) 

y-\\y-za 

= l + c- \\(p(t,t ,x) - z \\ . 
Now consider a linear system of differential equations 



dx 
~d~b 



A (t) x, 



(5) 



where A (t) is matrix of order n x n with continuous components, t G [to, T]. In 
such a case V (x, t) = 1 — c + \\X (tp, t) x — zq\\. Here X (t, to) is a fundamental 
matrix of solutions of the system (||) normalized with respect to to. The problem 
is to define the maximal radius c* of the boll K c (zo) such that th e tri vial solution 
of the system (||) is {-Kc* (^o) , ^t, to, T}-stable. From theorem 3J follows that 
we choose parameter c so that \lj (t) x\ < 1, s = 1, 2, . . . , N as soon as V (x, t) < 
1. Maximizing Lagrange function 

L (x, A) = ijx + A (\\X (t a ,t) x - zof - c 2 ) , s £ {1, 2, . . . N} , 
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t e [t ,T] yield 

l-\lJ(t)X(t,t ) zq\ 



mm mm 



te[t ,T]s=l,2,...N ,JlT ( t ) JJ i g (fi 



(6) 



where the matrix H (t) = X (t, to) X T (t, to) satisfies the Lyapunov's differential 
matrix equation 



dt 



= A(t)H(t)+H{t)A I (t),te[t ,T] (7) 



with initial condition H (to) — I, I is unit matrix of order n. 
Let Q be a positive definite symmetric matrix, 



K. 



(z a ) = |x € R n : (x - z ) T Q(x- z ) < c\ 



Consider the set of initial conditions in the form Go = K c (zq). It is necessary 
to obtain the maximal c* from such c > that the trivial solution of the system 

(||) is {K c (z ) , r t , t , T}-stable. If we choose p (a, b) = \J (a — b) T Q (a — b) = 
\\a — b\\g then V (x, t) = 1 — c + \\X (to, t) x — z \\q. Proceeding as it stated 
above we obtain 

i-\iT{t)x(t,t )z \ 

mm mm 



te[t,T] s =i,2,...JV ^ij ( t )H(t) l s (t) 
Here H (t) = X (t, t ) Q- X X T (t, t ) satisfies (@), H (t ) = Q~ x . 
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